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Abstract 

This paper is concerned with the study of a circular random distribution called 
geodesic Normal distribution recently proposed for general manifolds. This 
distribution, parameterized by two real numbers associated to some specific 
location and dispersion concepts, looks like a standard Gaussian on the real 
line except that the support of this variable is [0, 27r) and that the Euclidean 
distance is replaced by the geodesic distance on the circle. Some properties are 
studied and comparisons with the von Mises distribution in terms of intrinsic 
and extrinsic means and variances are provided. Finally, the problem of esti- 
mating the parameters through the maximum likelihood method is investigated 
and illustrated with some simulations. 

Circular statistics deal with random variables taking values on hyperspheres and can be in- 
cluded in the broader field of directional statistics. Applications of circular statistics are numerous 
and can be found, for example, in fields such as climatology (wind direction data |MJOO| ). biol- 
ogy (pigeons homing performances IWatSSI ) or earth science (earthquake locations occurence and 
other data types, see |MJOO| for examples) among others. 

A circular distribution is a probability distribution function (pdf) which mass is concentrated 
on the circumference of a unit circle. The support of a random variable 9 representing an angle 
measured in radians may be taken to [0, 27r) or [— tt, tt). We will focus here on continuous circular 
distributions, that is on absolutely continuous (w.r.t. the Lebesgue measure on the circumference) 
distributions. A pdf of a circular random variable has to fulfill the following axioms 

(i) .f{0) > 0. 

(ii) j^-fie)de^i 

(iii) f{9) = f{6 + 2kTT) for any integer k (i.e. / is periodic). 
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Among many models of circular data, the von Mises distribution plays a central role (essentially 
due to the similarities shared with the Normal distribution on the real line). A circular random 
variable (for short r.v.) 9 is said to have a von Mises distribution, denoted by vM{ii, k), if it has 
the density function 

2ttIq{k) 

where fi G [0,27r) and k > are parameters and where Io{k) is the modified Bessel function 
of order 0. The aim of this paper is to review some properties of another circular distribution 
introduced by |Pen06| which also shares similarities with the Normal distribution on the real line. 
For some fi G [0, 27r) and for some parameter 7 > 0, a r.v. 6 is said to have a geodesic Normal 
distribution denoted in the following gN{^,j), if it has the density function 

where dai', ■) is the geodesic distance on the circle and where k(-j) is the normalizing constant 
defined by 

/27r 1271 ( /t" 

^(7) := y — (^Uv^) - *(-7r/7) = Y — ( 2 

where $ is the cumulative distribution function of a standard Gaussian random variable and where 
er/(-) is the error function. 

Let us underline that Pennec |Pen06) introduced the geodesic Normal distribution for general 
Riemannian manifolds. We focus here on a special (and simple) manifold, the circle, in order 
to highlight its basic properties and compare them with the most classical circular distribution, 
namely the uM(/x, k) distribution. That is, we present here a new study of the gN distribution 
in the framework of circular statistics and provide results in terms of estimation and asymptotic 
behaviour. 

One of the main conclusions of this paper may be summarized as follows. While the von Mises 
distribution has strong relations with the notion of extrinsic moments (that is with trigonometric 
moments), we will emphasize, in this paper, that the geodesic Normal distribution has strong 
relations with intrinsic moments that is with the Frechet mean (defined as the angle a minimizing 
the expectation of ^(^(q;,^)^) and the geodesic variance. vM and gN distributions definition are 
closely related to respectively extrinsic and intrinsic moments, and we present their similarities 
together with dissimilarities. 

After introducing the gN distribution in Section [TJ we present a brief review on intrinsic and 
extrinsic quantities that allow characterization of distributions on the circle in Section[2] In Section 
[3l we present extrinsic and intrinsic properties of vM and gN distributions. Then, in Section |4l 
we rapidly explain how to simulate gN random variables on the circle. Finally, in Section O we 
present the maximum likelihood estimators for the gN distributions and study their asymptotic 
behaviour. Numerical simulations illustrate the presented results. 
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1 Geodesic Normal distribution through the tangent space 



As introduced in |Pen06| the geodesic Normal distribution is defined for random variables taking 
values on Riemannian manifolds and is based on the "geodesic distance" concept and on the total 
intrinsic variance |iBP03[ |BP05[ [Pen06| . On a Riemannian manifold A^, one can define at each 
point X E A4 a scalar product (., .)^ in the tangent plane T^Ai attached to the manifold at x. 
On A4, among the possible smooth curves between two points x and y, the curve of minimum 
length is called a geodesic. The length of the curve is understood as integration of the norm 
of its instantaneous velocity along the path, and with the norm at position a; on taken as: 
II -11 = V ■)x- well-known that, given a point x E A4 and a vector if G T^M, there exists 
only one geodesic "f(t) with "f{t — 0) — x and with tangent vector if . 

Through the exponential map, each vector S T^Ai is associated to a point y £ A4 reached 
in unit time, i.e. 7(1) — y. Using the notation adopted in jPen06| . the vector defined in T^M. 
associated to the geodesic that starts from x at time i = and reaches y at time i = 1 is denoted 
xi). Thus, the exponential map (at point x) maps a vector of T^M to a point y, i.e. y — exp^(x^). 
Now the geodesic distance, denoted dG{x,y), between x A4 and y e is: 



dG{x,y) = \J {xfi,xfj)^ 



The Log map is the inverse map that associates to a point y G in the neighbourhood of a; a 
vector xfj e T^M., i.e. xfj — Logx{y). 

A random variable Y taking values in M.. with density function f{y;fj,,T) is said to have 
a geodesic Normal distribution with parameters fi G A4 and F, a {N, N) matrix, denoted by 
g7V(^,r),if: 

/(y;M,r) = fc^^exp ( ' ''^ 
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The parameter n is related to some specific location concept. Namely, |Pen06) has proved that 
H corresponds to the intrinsic or Frechet mean of the random variable Y (see Section [2] for more 
details). The normalizing constant is given by: 

fc = /^exp(-^^)dA^(,) 

where dA4{y) is the Riemannian measure (induced by the Riemannian metric). The matrix F is 
called the concentration matrix and is related to the covariance matrix of the vector fiY given for 
random variables on a Riemannian manifold by 



dM{y). 



Note that in the case where A^ = R'^, the manifold is fiat and the geodesic distance is nothing 
more than the Euclidian distance. In this case, we retrieve the classical definition of a Gaussian 
variable in R'' with ji and F corresponding respectively to the classical expectation and to the 
inverse of the covariance matrix S^^. 
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The case of the circle We now present, as done in |Pen06| . the case wliere Ai is the unit circle 
S^. The exponential chart is the angle 6* e] — tt;tt[. Note that this chart is "local" and is thus 
defined at a point on the manifold. This must be kept in mind especially for explicit calculation. 
Here, as the tangent plane to is the real line, 9 takes values on the segment between — tt and 
TT. Note that tt and — tt are omitted as they are in the cut locus of the "development point" (point 
on the manifold where the tangent plane is attached). dAi is simply d9 here. 

As stated before, the gN distribution on the circle has density function given for S {j-i ~ 
7T,n + Tr) by: 

where 7 is a nonnegative real number. Note that dclM, ^) is the arc length between /i and 9. The 
normalization k{'y) is: 

Hi) = l"^\-i'^oi,,sr^9 = ^erf (tt^I) , 

with the development made around /i. 

In order to consider a gN distribution as a classical circular distribution (defined by axioms (i)- 
(iii) in the introduction), one must extend the support of this distribution from (/i — tt, /i + 7r) to M. 
The way to achieve this is to make the geodesic distance periodic. Let us consider the distance da 
for an angle a G M \ + /ctt, fc e Z} defined by da{^, a) = dGi^J■^ 5) with a = (a — /i + tt) (mod 2tt). 
Let / the density / where da is replaced by da- This new density defines a circular disribution 
satisfying axioms (i)-(iii). In particular, 

/ f{9,fi,-f)d9= [ A:-i(7)e-*'^~'5(^^'')'(i6l = fc-i(7)e-^''<=('^''')'d6' = 1. 

Jo Jo J n-Tx 

For the sake of simplicity, da will be understood as the distance da in the rest of the paper. And 
therefore, the density of a gN distribution is considered periodic, defined in R\{/i + fc7r, fc G Z} and 
with values on M. In the case of a vM distribution (and actually for most of circular distributions) 
no such considerations are needed; the periodic nature being included through the cos(-) function. 

2 Classical measures of location and dispersion for circular 
random variables 

We briefly present the concepts of extrinsic and intrinsic moments for random variables on the 
circle. While the former are well-known in circular statistics [e.g. |MJ00| 1. the later, based on 
the geodesic distance on the circle, are less used in this domain. They have been introduced and 
commonly used when dealing with genereal Riemannian manifolds {e.g. |Kar77| [Zie77|, lBP03| and 
the numerous references therein). 
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Extrinsic moments 



In circular statistics, it is well established that trigonometric moments give access to measures 
of mean direction and circular variances. Considering a cicular random variable 6, its p — th 
order trigonometric moment is defined as ipp = E[e'^P^] = ap + i/3p where Up = E [cos p9] and 
/3p = E[sinp9]. These later quantities are extrinsic by definition. The first order trigonometric 
moment is thus ipi = i?[e*^] = pe^'^ where p is called the mean resultant length (0 < p < 1) and 
/i^ is the mean direction. In the following, we refer to p,^ as the extrinsic mean. The extrinsic 
variance cr|, is indeed the circular variance defined as a'^ = 1 — p. In the sequel, extrinsic moments 
will be used in place of trigonometric moments, keeping in mind that they are the same quantities. 
For more details on trigonometric moments see |MJOO| PASOI) . 

Intrinsic moments 

Another way to consider moments for distributions of random variables on the circle is to use the 
fact that is a Riemannian manifold and thus the geodesic distance can be used to define intrinsic 
moments. Given a random variable 9 with values on , we define by p^ = argmin^g^i {p, 0)^] , 
the intrinsic mean set (a particular case of a Frechet mean set), where we recall that dci-, ■) is the 
geodesic distance on the circle, i.e. the arc length. When the intrinsic mean set is reduced to a 
single element, p^ is simply called the intrinsic mean. The intrinsic variance is then uniquely 
defined hy aj = E [^^(/i^, 9)"^] where p^ is the intrinsic mean (set) defined above. For more details 
and a thourough study of intrinsic statistics for random variables on Riemannian manifolds, see 
|Pen06) or |BP03[[BP05| . Other concepts of intrinsic variance (e.g. variances obtained by residuals 
and by projection) exist (see |HHM10) for thourough description). In this paper, we only focus 
on the intrinsic mean and (total) variance which are sufficient to basically compare gN and vM 
distributions. 

3 Basic properties of the gN{ii,^) and vM{fi,K) distributions 

Symmetry property 

First, let us say that like the vAI distribution, the gN distribution has a mode for 9 = p, and 
a symmetry around 9 — p. The vM distribution has an anti-mode for 9 ~ p ± tt. For a gN 
distribution, the density is not defined at these points. However, the shared behaviour is the 
decreasing of both densities on each interval (p — tt, p] and [p, p + tt). 

Extrinsic and intrinsic means and variances 

Table [1] summarizes extrinsic and intrinsic means and variances for both distributions of interest. 
Let us make some comments. In [KSOSj , the authors follow the works of Le ( |Le98| ILeOlj ) and 
give very simple conditions on the density of a circular random variable that ensure the existence 
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and unicity of the intrinsic mean. It is left to the reader that applying Theorem 1 of |KS08| allows 
us to assert that the intrinsic mean of a gN{fi, 7) or vM{ji^ k) is /i. Furthermore, the computation 
of for a gN distribution (resp. a'^ for a vM distribution) can be found in |Pen06| (resp. 6.17. 
|MJOO| ). The intrinsic variance trj for a vM is quite obvious and is omitted. It remains to explain 
how we obtain the extrinsic moments for a gN distribution. Both are indeed derived through the 
pth trigonometric moment of a gN distribution reported in the following proposition. 

Proposition 1. The p-th trigonometric moment (p G N* ) of a gN(fj,,"f) distribution, denoted by 
(fip and defined by ipp :— E [e*P^] is given by 

^2 Reierf U^-i-^ 



where erf is the error function defined for any complex number by erf{z) = /-^oo 
Proof. Let p > 1, 

Since for e (/i — tt, + tt), dciO, /i) = \0 — 
E[e'P^] = fc-i(7) / e'P^^+^'^e-i^'d9 



— IT 



— 7r 



27r 1 / / /t" ■ P \ r f pi P 



e 2^ 



72V V V 2 V27/ V V 2 ^/2^ 

Re { erf ^vr-^/y — i 



2 727 



since for any complex number z, erf{z) — erf{—z) = erf(z) — erf{—z) = erf{z) + erf{z) = 
2Re{erf{z)). U 

Distribution fx^ jj,^ aj a"^ 

g^(M,7) l(l-2.fc-^(7)e^^) ^ - 



Table 1: Summary of extrinsic and intrinsic means and variances for the von Mises and geodesic 
Normal distributions. The function Ip{ ) denotes the modified Bessel function of the first kind of 
order p. 
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Figure [3] shows the evolutions of the extrinsic and intrinsic variances in terms of the concen- 
tration parameter k for the vM and 7 for the gN. It is interesting to notice that the von Mises 

2 

distribution and the geodesic Normal distributions have intrinsic variance equal to ^ when 7 or 
K equals zero, corresponding to the variance of the uniform distribution on the circle. Note also 
that intrinsic and extrinsic variances tend to zero as 7 or k tend to infinity. 




Concentration parametEr 



Figure 1: Evolutions of extrinsic and intrinsic variances in terms of the concentration parameter, 
7 for the gN distribution and k for the vM one. 

Entropy property 

gN and vAI have both a peculiar position respectively amongst distributions defined on Rieman- 
nian manifolds and circular statistics distributions: both maximize a certain definition of the 
entropy. As explained in [JASOlj (characterization due to Mardia [Mar72| l. the circular distribu- 
tion that maximizes the entropy defined with respect to the angular random variable 6 (see |iMJOO| 
for exact definition), subject to the constraint that the first trigonometric moment is fixed, i.e. for 
IJ,^ and a'^ fixed, is a vM(/i^.K), where a'^ = 1 — /i(k)//o(k). In a similar way, as demonstrated in 
|Pen06) ■ the distribution defined using the geodesic distance on which maximizes the entropy, 
when it is defined in the tangent plane, and subject to the constraints that /i^ and aj are fixed, is 
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the gN{fi^ One can thus conclude that vM and gN distributions play "similar" roles in that 
they maximize the entropy with respect to either extrinsic or intrinsic moments of the distribution. 

Linear approximation of /i0 

Recall that the random variable /i0 represents the algebraic measure of the vector fj.9. The support 
of Iji9 is {—TT, tt). Its cdf is given for t e (— tt, tt) by 

This reduces to F—K(t) — } { — 7 — ^^-4-. In other words, ud is nothing else than a truncated 
Gaussian random variable with support (— tt, tt) with mean and scale parameter \l ^J^, that is 

= Z| |Z| < TT, where Z-7V(0,1/V7)- (1) 

For large concentration parameters 7, the geodesic Normal distribution can be "approximated" 
by a linear normal distribution in the following sense. 

Proposition 2. Let 6 ^ gN{fi,j), then as 7 — > +00, .^7 fi9 -4 J\f{0, 1). 

Proof. Let us denote by h{-) the moment generating function of the random variable ^/7 ^9, 
then for i e M 

h{t) = k-\j) [ e^'^e-i^'dO 

J —TT 

= k-\^)^ e-i^'-'^'d0 

^ $ (7r77 ~ ^) - $ (-TT^ - t) 
$ {tT^) - $ (-TT^) 

Therefore, as 7 +00, for fixed t, h{t) converges towards e* which is the moment generating 
function of a standard Gaussian random variable. H 

Such a result also holds for the von Mises distribution: as k ^ +00, ^/k,{6 — ^) A- J\f{0, 1), see 
e.g. Proposition 2.2 in jJASOl) . 

4 Simulation of a geodesic Normal distribution and examples 

The generation of a gN distribution with support on (0, 27r) is extremely simple following ([T]). It 
consists in two steps. 

1. Generate Z\ \Z\ < tt, where Z ^ J\f{0, 1/77) 

2. Set 6* = ^ + Z(mod 2tt). 

Figure [2] presents some examples for different values of location parameter /x and concentration 
parameter 7. 
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Figure 2: Example of n = 100 gN{fj,,^) realizations for different parameters. The red line indi- 
cates the value of /i. These different plots have been produced using the R packages circular, 
CircStats maintained by C. Agostinelli, and rej^ted to the book [JAS01| . 



5 MEiximum Likelihood Estimation 



5.1 Preliminary and notation 

Let us consider now the identification problem of estimating the parameters of a gN distribution 
from the n observations 0i,. . . ,9n G (0, 27r). In this section, we will denote by fj,* and j* the 
unknown parameters to estimate. We assume that fi* G [0, 27r) and 7* > 0. Also, we propose to 
denote by fl^ and /I^ the empirical intrinsic and extrinsic means defined by 

1 " 

Ji^ := argmin^g5i-^dG(/U,6'i)2. (2) 
" i=i 

:= Arg(^i), with ^1 := — cos((^j) + i— sin((^j). (3) 



n n 
3 3 



Obtained through the minimization of an empirical function, ju^ is not necessarily reduced to a 
single element. The natural intrinsic and extrinsic variances are then denoted by a| and d% and 
uniquely given by 

1 " 

d'i = -^dcill^ ,eif and ct| = 1-|^i|. (4) 



" • 1 

z=l 



In the following, we will need the following Lemma and notation. 

Lemma 3. For a random variable ~ gN{^*,^), let ^(/Lt, 7) := E[dG{lJ',Oj)^] for ji G (/x* 
TT, /X* + tt) and 5 = ji* — ji, then 

^(m,7) =ff('5)l[o,x)('^) +5(-'5)l(-.,o]('^), 

where 

/TT pTT 
{S + aff{a)da + An j (27r -{5 + a))f{a)da 

and f{a.) = k~^{j)e~i°'^ . 

Proof. Let us fix G (/x* — tt, /x*] , then 

{fifi* + fj*a when a € (fi* — tt, fi* + tt -~ (/i* — /i)) 

IJH* + iJ*a — 27r when a £ {/i* + n — {fj.* — /x), n* + tt) 

Denoting 6 := jj.* — ji, this expansion allows us to derive 

y(/x,7) = E[dGi^i,0^)'] 



{6 + aff{a)da+ / {2Tr - 6 - af f {a)da 

-TV J-K — 5 

/TT pTT 
{S + aff{a)da + 4n (27r -{6 + a))f{a)da. 
-TT JtT — 6 
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Now, let II G [fi* , fi* + tt) , then 
/IS = 



/i/i* + j-i*a when a G {p* — vr + (/t — fJ.*), fi* + tt) 

/i/i* + i-i*a + 27r when a G (/i* — tt, /.t* — tt + (/i — //*)), 



which leads to 



-TT — (5 

F(Ai,7) = / {6 + a)^f{a)da+ I {2tt + 6 + af f{a)da 

-7r — 5 J —TT 

/ — 7T — S 
{2Tr + {S + a)) f{a)da 
-TT 

(S + aff{a)da + 4tt {2tt + {6 ~ a))f{a)da 

T JtT + S 

■ 

Obviously V{^*,'y) corresponds to the intrinsic variance of 9^. Then, from Table [1] this 
function does not depend on /i* and will therefore be simplified to Vlj). As used in Theorem 21 
let us recall here the expression of the later quantity. 



y(7) -F(/.*,7) := - (^l-2^fc-i(7)e- — j with ^(7) = J— er/ ( ttJ^ ) (5) 

5.2 Maximum Likelihood Estimate 

The log-likelihood expressed for n i.i.d. gN distributions is given by : 



£(m,7) = -nlog(fc(7)) ~ ^Y.'^GilJiAf ■ 

1=1 

Let = {(/x,7) : /X G [0, 27r) \ {(/i* ± tt) mo(i(27r)},7 > 0} and assume that the true parameter 
(/i'^,7*) belongs to the interior of il. The MLE estimates and asymptotic results are given by the 
following result. 

Theorem 4. 

(i) The MLE estimate of fi* corresponds to the intrinsic sample mean set, that 

The MLE estimate of j* is uniquely given byj^^^^ — V^^{aj), where V{-) is the function defined 

by (0) and where aj is the intrinic sample variance. 

(m) As n — > +00, (/i*^^^, 7*^^^) is a strongly consistent estimate of {ji*,^*). 
{Hi) As n — > +00, the MLE estimates satisfy the following central limit theorem 

{^'^^ - M^ 7''^^ - I'^f ^ (0, j-i(7*)) , 

where J {'~f*) is the Fisher information matrix given by J {'J*) = | | with Ji{'y*) := 

V J2{j*) / 



7^ (1 - 27rfc-i(7*)e-^^') and 72(7*) := 
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We emphasize that we do not derive analytic formulas for but we prove its uniqueness by 

proving that the function V{-) is a strictly decreasing function (which is illustrated by Figure [3]). 
From a practical point of view, the computation of 7*^^^ (as well as the one of /I^^^^) has been 
derived using a simple optimization algorithm. 

As Mardia and Jupp did for the von Mises distribution ( |MJOO| . Section 5.3 p. 86), /i^^^^ is 
regarded as unwrapped onto the line for the asymptotic normality result. 

As it is for a Gaussian distribution on the real line or for a vM distribution on the circle, the 
Fisher information matrix of a gN distribution does not depend on the true location parameter 
IJ* and the two estimates of /i* and 7* are asymptotically independent. Let us also note that the 
geodesic moment estimates, that is the estimates of fi* and 7* based on the first two geodesic 
moments equations and V{j) exactly fit to the maximum likelihood estimates. Here is again 
another analogy with the vM distribution, since the MLE of a vM distribution correspond to the 
estimates of /i and k based on the extrinsic moments (see [MJOOj for further details). 
Proof. (i) Since the minimum of X]r=i '^g(M: ^i)^ defines the intrinsic sample mean set, the 
MLE of fi correponds to the intrinsic sample mean set of fi* . Now, the partial derivative of £ with 
respect to 7 is given by 



Let us note that 



^'^^^ = ^2 J / ' de = -^E[dG{^J^\9,Y] = -^V{^). (6) 

Replacing ^ by its MLE estimate and taking the derivative of £ w.r.t. 7 equal to zero implies that 
the MLE estimate of 7 is defined by the following equation: 



l^(7Mi^) = iVdG(/i^"^ 

n ^-^ 



1=1 



The proof is ended by showing that V{-) is a striclty decreasing function on R+. Similarly to ([S]), 
we notice that ^"(7) = ^E[dG{^J* ,9^)% Now, 

'1 



= -2 \^-E[dG{^i\e^t] - E[dG{^i\ o^y/i 

= -iFar[o^G(/i^^)'] <0. (7) 

[a) As fjL* corresponds to the unique intrinsic mean for a gN distribution, the strong consistency of 
the intrinsic sample mean set is derived from e.g. Theorem 2.3 of |BP03| . Now, for the consistency 
of 7^^^, let us introduce the variable M„ := X]r=i dG{lJ'* , di)"^. From the LLN, M„ converges 
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almost surely towards i?[c?G(/i*, 0^)^] — ^^(7*). Moreover, we may prove that A/„ — ct^ tends to 
zero (almost surely). The later following from 



\M^-d\\ 



n 



in 



< —y^\dG{e^,^l 



MLE\ 



47r , 

< — 



n 



MLE *\ 
I M ) 



1=1 



Combining the previous convergences leads to the almost sure convergence of af to ^^(7*) and to 
the result since V{-) is continuous and invertible on W^. 

{in) Standard theory of maximumn likelihood estimators (Theorem 5.1 p. 463 of |LC98[ ') shows 
that asymptotic normality result holds. The verification of the assumptions (A-D) of |LC98) . 
p. 462-463 are omitted; we just focus, here, on the computation of the Fisher information matrix. 
The antidiagonal term is given by 



J12 :- -E 



d_ 



dGifJ;dyy 



IdV 
2 9^ 



(a**,7*)=0. 



since fi* corresponds to the intrinsic mean and thus minimizes the geodesic variance. The asymp- 
totic variance of y^7^^^^ is given by the inverse of 



^2(7*) = 



fc"(7*)fc(7*) - fc'(7*)2 



*\2 



Recall that from ([7]), this constant is positive. Now, the last term to compute is the asymptotic 
variance of ^Jnfl^'^^^ given by the inverse of 



^1(7*) 



Ie 



g2 

-^dGiii,9^y 



2 



From Leinma[3l cr|(/i, 7) is a function of (5 = /i* — /i. Without loss of generality, assume ^ > (the 
other case leads to the same conclusion), then the function g (in Lemma [S]) is twice continuous 
differentiable on [0, tt) and g"{S) = 2 — 47r/(7r — S). Setting (5 = in the last equation leads to the 
stated result. ■ 



5.3 Simulation study 

We have investigated the efhciency of the maximum likelihood estimates in a simulation study. A 
part of the results are presented in Table [21 As expected, the empirical MSE of both estimates 
of the parameters /i* and 7* converge towards zero as the sample size grows. We also notice that 
it's more complicated to estimate the intrinsic mean when the concentration parameter is low. 
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Unlike this, the concentration parameter is better estimated for low values of 7*. These facts are 
confirmed by Figure [3] which shows the constants of the asymptotic variances (for both estimates), 
i.e. 1/Ji(7*) and 1/J2(7*), in terms of 7*. Figure U] illustrates the central limit theorem satisfied 
by the MLE estimates. 



Location parameter /i* 



Simulation 
Parameters 



Sample size 

n=m n = 20 71 = 50 n = 100 n = 500 







5 


3.3668 


1.6515 


0.2218 


0.0313 


0.0048 


fl* = 




1 


0.1088 


0.0537 


0.0206 


0.0103 


0.0020 




5ZL ^* - 
4 ' ' 


5 


0.0199 


0.0102 


0.0040 


0.0019 


0.0004 




ZZL ^* - 


10 


0.0099 


0.0050 


0.0020 


0.0010 


0.0002 



Concentration parameter 7* 



Simulation 








Sample size 






Parameters 




n = 10 


n = 20 


n = 50 


n = 100 


n = 500 


M* = f,7* = 


5 


0.2893 


0.0748 


0.0255 


0.0141 


0.0065 




1 


1.0559 


0.2268 


0.0565 


0.0252 


0.0046 




5 


22.8649 


5.4092 


1.3365 


0.5757 


0.1037 




10 


88.5093 


20.1554 


5.3607 


2.2723 


0.4179 



Table 2: Empirical Mean Squared Error (MSE) of MLE estimates of the location parameter 
fi* (top) and the concentration parameter 7* (bottom) based on m = 5000 replications of gN 
distributions for differents choices of parameters and different sample sizes. 
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Figure 3: Plots of the constants of the asymptotic variances of /i^^^^ (left) and j'^'^^^ (right), i.e. 
the constants 1/Ji(7*) and 1/J2(7*) given in Proposition lU in terms of 7*. 
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Figure 4: Histograms of (y/n(jlf^'"^ — fi*)) . based on m = 5000 replications of a 

^ J 'J — J. , . . . , DUUU 

(7iV(37r/4, 1) distribution for different sample sizes. The curve corresponds to the density of a 
Gaussian random variable with mean zero and variance 1/Ji(l) given in Proposition 2] 
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